Na2Vs07, a frustrated nanotubular system with spin-1/2 diamond rings 
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Following the recent discussion on the puzzling nature of the interactions in the nanotubular 
system Na2V307, we present a detailed ab-initio microscopic analysis of its electronic and magnetic 
properties. By means of a non-trivial downfolding study we propose an effective model in terms 
of tubes of nine-site rings with the geometry of a spin-diamond necklace with frustrated inter-ring 
interactions. We show that this model provides a quantitative account of the observed magnetic 
behavior. 



Introduction- Low-dimensional quantum spin systems 
with chain, ladder or planar geometries have attracted 
much attention in the last yearsi due to their uncon- 
ventional magnetic properties. Recently substantial ef- 
fort has been devoted to the search for new materials 
with exotic topologies and novel properties. For instance, 
the group of Millet et al£ has succeeded in synthesizing 
the first transition-metal-oxide-based nanotubular sys- 
tem Na2 V3O7. The discovery of Na2 V3O7 has introduced 
a new class of novel geometric structures in the field of 
magnetic nanostructures. 

The V 4+ C>5 pyramids in Na2V307 share edges and cor- 
ners to form a nanotubular structure - a geometry first 
discovered for carbon- with Na atoms located inside and 
around each individual tube (see FIG. QJ. The complex 
geometry of the compound is expected to provide non- 
trivial paths for the exchange interaction and hence to ex- 
hibit a non-trivial magnetic behavior. Initially this com- 
pound was predicted to be an example of a S = 1/2 nine- 
leg (or three-leg ladder) system with periodic boundary 
conditions along the rung direction^. Using the extended 
Hiickel method, Whangbo and Koc 3 conjectured, on the 
other hand, that the tubes could be described by six mu- 
tually intersecting helical spin chains which should show 
a gap in the spin excitation spectra. 

Unlike antifcrromagnctic (AF) spin ladders with even 
number of legs, the behavior of odd- leg ladders is strongly 
influenced by topology. An odd-leg ladder with open 
boundary conditions in the rung direction behaves effec- 
tively as a S — 1/2 AF Heisenberg chain and therefore 
shows no gap in the spin-excitation spectra. This situa- 
tion changes dramatically when the boundary condition 
along the rung direction is periodic, as is the case for 
the spin tubes in Na2Vs07. The introduction of peri- 
odic boundary conditions brings in an additional degree 
of freedom in the problem, namely the chirality^*£*£*£ and 
the ground state is then four-fold degenerate with a gap 
to the first excited state. However, the magnetic proper- 
ties of Na2Va07 as reported in the literature neither fit 
the properties of odd-leg spin tubes nor to that of helical 
spin chains as proposed by Whangbo and Koo^ since no 
appreciable spin gap could be detected in the susceptibil- 
ity measurements^. A way out of this puzzling situation 
has been suggested by Liischer et. al& by considering 



frustrated inter-ring couplings within the framework of a 
three-leg spin tube. While this surely provides an inter- 
esting model, its relevance for Na2Va07 remains yet to 
be examined. 

Eventhough this compound has generated interest and 
controversy, no first-principles microscopic study has 
been reported so far. Our work in that respect is 
the first microscopic study carried out for Na2V307. 
We have analyzed the ab-initio density functional the- 
ory (DFT) results in terms of the newly developed N- 
th order muffin-tin-orbital (NMTO) based downfolding 
technique^ to provide a microscopically-derived spin 
model for Na2V307. Our results show that the appro- 
priate description of the system is that of tubes of nine- 
site rings with partially frustrated, next-nearest neighbor 
intra-ring interactions and frustrated inter-ring interac- 
tions of weaker strength. Validity of this model is pro- 
vided by the good agreement of our calculated suscepti- 
bility data with the measured data. 

Ab-initio Study- Na2V307 crystallizes^ in the trigonal 
space group, P31c. The unit cell contains six formula 
units and the lattice constants are a — 10.89 A and c = 
9.54 A. There are three inequivalent V sites -VI, V2, V3, 
seven inequivalent O sites - 01-07 and four inequivalent 
Na sites, Nal- Na4. Each V is surrounded by five oxygen 
atoms with one characteristic short V-0 bond giving rise 
to distorted square pyramids which are connected to each 
other to form nanotubes with internal diameter of about 
5 A (see FIG. EDO). The 

arrangement of VO5 pyramids 
in a single tube is better viewed in the unfolded idealized 
representation of the nanotube as shown in FIG. db). 
The 01, 02 and 03 atoms, not marked in the figure, 
are positioned at the apex of VO5 pyramids centered by 
VI, V2 and V3 atoms respectively, all of them point- 
ing out of the tube. The edge sharing VO5 pyramids 
form nine-member rings out of the basic unit V1-V2-V3, 
in alternating sequence of (V2-V1-V3)-(V2-V1-V3)-(V2- 
V1-V3) and (V3-V1-V2)-(V3-V1-V2)-(V3-V1-V2) within 
slices (A) and (B). Rings in slices (A) and (B) are con- 
nected to each other by the corner sharing oxygens 05 
and 06 to form a nanotube oriented along the crystallo- 
graphic c-direction. 

The DFT band-structure in the generalized gradi- 
ent approximation (GGA)ii, calculated using the lin- 



2 





FIG. 1: (Color on-line) Crystal structure of Na2Vs07- (a) 
Projection on to (001) plane. V, Na and O atoms are denoted 
by black, blue and red balls respectively, (b) Unfolded view 
of the edge and corner-sharing VO5 pyramids within a tube. 



earized muffin tin orbital (LMTO) method based on 
the Stuttgart TBLMTO-47 code±2 and the linear aug- 
mented plane wave (LAPW)i£ method, consists of pri- 
marily oxygen— p and V— d derived bands in the region 
of interest^. The oxygen— p dominated bands are sep- 
arated by an energy gap of about 3 eV from the V— d 
dominated bands . The V— d dominated bands span an 
energy range of about 4 eV, starting from about -0.7 eV 
to 3.5 eV, with the zero of energy set at the GGA Fermi 
energy. The pyramidal co-ordination of oxygen atoms 
surrounding the V atom, sets the V—xy d orbital to be 
the lowest energy state within the V— d manifold (in the 
local reference frame, with the z-axis pointing along the 
shortest V-0 bond and the x-axis pointing along the next 
shorter V-0 bond). The V—xy bands therefore appear 
as essentially non-bonding set of bands extending from 
about — 0.7eV to 0.4 eV around the Fermi level, with 
small mixing of oxygen character. The crystal-field split 
V— yz, xz, 3z 2 — 1 and x 2 — y 2 complexes appear in the 
energy spectrum above the V—xy complex in order of in- 
creasing energy. Na-derived states lie farther high up in 
energy with practically no mixing to bands close to the 
Fermi energy^. 

Starting from such a density- functional input, it is a 
non-trivial task to build up a low-energy model Hamil- 
tonian relevant for the system. For the sake of unique- 
ness, however, it is essential for such model Hamiltoni- 
ans to be derived in a first principles manner containing 
the essential chemistry of the material. In recent years 
the N-th-order-muffin-tin-orbital-method-based™ down- 
folding technique has been successful in achieving this 
goalie. The method relies on designing energy-selective 
Wannier-like effective orbitals by integrating out degrees 
of freedom that are not relevant - a method called down- 
folding. The few-orbital Hamiltonian is then constructed 
in the basis of these Wannier-like effective orbitals. Since 
the degrees of freedom are integrated out, rather than 
being simply ignored, the method naturally takes into 




FIG. 2: Hopping parameters for Na2Va07 for V atoms be- 
longing to slices (A) and (B) in FIG. 1(b). For identification 
of V atoms belonging to slice (A) refer to FIG. 1(b). 



TABLE I: Inter and intra-ring hopping integrals in eV. The 
next important hoppings are about 0.001 eV. 



Intra-ring 


Inter-ring 


tl tl t-$ 


t'i t' 2 


tl± t2± 


*1_L *2X 


-0.18 -0.15 -0.14 


-0.13 -0.14 


-0.03 -0.02 


-0.02 -0.03 



account the renormalization effect coming from orbitals 
not considered explicitly. In particular, in the present 
case, we integrate out all the degrees of freedom other 
than V—xy orbitals. The effective V—xy muffin tin 
orbitals (MTO's) generated in the process contains in 
its tail the integrated out O— p and remaining V— d or- 
bitals, the weight being proportional to their mixing to 
V—xy derived bands. Fourier transform in the down- 
folded V—xy basis gives the tight-binding Hamiltonian, 

Htb — — 52<i j>tij( c j c i + c l c j) m terms of dominant 
V-V effective hopping integrals, iy, where i and j denote 
a pair of V 4+ ions. 

The nearest neighbor (n.n.) V-V interactions within 
the V1-V2-V3 basic unit (marked by dark green in 
FIG. Gib)) proceed via the edge-sharing oxygens while 
the intra-ring next nearest neighbor (n.n.n.) interac- 
tions proceed via corner-sharing oxygensit Our first- 
principles-derived hopping integrals show that the edge- 
sharing n.n. V-V interactions within the V1-V2-V3 basic 
unit, ti, t 2 and t 3 (see FIG. 2 and TABLEQ}, are of mag- 
nitude ranging from -0.14 eV to -0.18 eV. The n.n.n cor- 
ner sharing intra-ring interactions t\ and t' 2 are nearly 
equally strong. The inter-ring V-V couplings ti±, t' i± 
z = l,2 are an order of magnitude weaker than the intra- 
ring couplings. Nevertheless, we notice that the n.n.n 
inter-ring interactions (t 1± and t 2 j_) are equally strong 
as the n.n. inter-ring interactions (ti± and i2±). This 
induces inter- ring frustration which, as we will see later, 
could be of fundamental importance for the description 
of the magnetic behavior of the system at low tempera- 
tures. 

This result is very different from the extended Hiickel 
molecular-orbital-based result of Whangbo and K00& 
which predicts that coupling via corner-sharing pyramids 
is much larger than that via edge-sharing for Na2Vs07. 
While this conclusion is in general true for vanadate sys- 
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FIG. 3: (Color on-line) Various intra-ring and inter-ring over- 
laps of the downfolded V—xy MTO-s at two VI and V2 sites. 
Plotted are the orbital shapes(constant-amplitude surfaces) 
with the lobes of opposite signs labeled by black(magenta) 
and white(cyan) at site VI (V2). (a) Edge-sharing intra-ring, 
nearest-neighbor VI- V2 overlap (ti) (b) Corner-sharing intra- 
ring, second neighbor VI- V2 overlap (t'i) (c) Inter-ring VI- V2 
overlap (tix)- 



terns like NaV 2 5 , LiV 2 5 and CsV 2 03±& it fails in the 
case of Na 2 V307. Our result shows that the peculiar ge- 
ometry of the VO5 coordination with marked deviation 
of the V-O-V bond angles from 180° and highly 11011- 
coplanar nature of V-O-V bonds, gives rise to coupling 
via edge-sharing nearly equally strong as that via corner- 
sharing as long as one is confined to a single ring. The 
magnitude of intra-ring, edge-sharing couplings are com- 
parable to the largest edge-sharing hopping parameters 
in other vanadate systems like CsV 2 Os (-0.12 eV)i^ or 
LiV 2 05 (-0.18 eV)2£ while the magnitude of intra-ring, 
corner-sharing couplings is much smaller than the largest 
corner-shared V-O-V in LiV 2 05 or NaV 2 0s which are 
about -0.3 e\S2*2i. The corner-sharing coupling between 
two adjacent rings is substantially diminished due to fur- 
ther mis-alignment of V—xy orbitals belonging to two 
different rings which mixes the tt character of the bond 
with that of S character. In FIG. we show the var- 
ious downfolded V—xy MTO-s and their overlap giving 
us an idea of the influence of the distorted geometry on 
the relative orientations of the effective V—xy orbitals, 
interaction paths and the magnitude of overlaps. 

Our microscopically derived hopping integrals show 
that Na 2 Va07 can be described as formed by tubes con- 
sisting of frustrated weakly coupled nine-site rings with 
partial intra-ring frustration. The partial frustration is 
due to the absence of the second neighbor VI- V3 interac- 
tion in the ring which is neither corner-sharing nor edge- 
sharing, rather it is decoupled by the intervening V2(0^) 
pyramid. 

Susceptibility - The exchange integral, J, can be ex- 
pressed in general as a sum of antiferromagnetic and fer- 
romagnetic contributions 22 J = J AF + J FM . In the 
strongly correlated limit, typically valid for transition 
metal oxides, the antiferromagnetic contributions are re- 
lated to the hopping integrals, U, by using a perturbative 
approach as J AF — 4tf / (U —Vmter) , where U corresponds 
to hopping via various V-O-V superexchange paths and 
U and Vi n ter are the on-site and intersite Coulomb in- 
teractions respectively. In the absence of a generally 



accepted satisfactory way of direct computation of ex- 
change integrals, such an approximate method is a good 
starting point for estimates of the exchange couplings as 
well as the relative strengths among the various exchange 
integrals 2 ^. 

In view of our ab initio results, the underlying spin-1/2 
Hamiltonian for Na 2 Va07 can be written as 

9 

H = Ji ^ ' + a t Si ■ Si+2) , (1) 
1=1 

where we neglect the small differences in between the 
three n.n. and two n.n.n intra-ring interactions and to a 
first approximation also neglect the small inter-ring cou- 
plings, cti — J\j J\ and we impose periodic boundary 
conditions S^-H = Si (L=9). In this model, the n.n.n 
coupling is inhomogeneously distributed in the sense that 
on = for i = 1, 4, 7 and a* 7^ (=«') for i = 2, 3, 5, 6, 
8, 9 (see left illustration in FIG. 0}. For a further check 
of the goodness of this model we considered in addition 
the fully frustrated model with J\ = J 2 for all i (a, = a) 
(see right illustration in FIG.0J. 
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FIG. 4: Temperature dependence of the inverse magnetic sus- 
ceptibility in units of mol/emu obtained from exact diagonal- 
ization for the two frustrated models of the inset (see Eq. Q) 
compared with the experimental data— (filled dots). 

In FIG. 0] we present a comparison of the experi- 
mentally observed inverse magnetic susceptibility- with 
that obtained from exact diagonalization of the above- 
mentioned two kinds of frustrated models, varying the 
parameters J\ as well as a or a'. 

We observe that only the partially frustrated model 
(ct^= 0), which is the model predicted from the ab initio 
calculations with J\ and J 2 antiferromagnetically signed 
is able to reproduce the experimental data over the whole 
range of temperature. The fully frustrated model consis- 
tently shows, on the other hand, an upturn of the inverse 
susceptibility at lower temperatures not observed exper- 
imentally. We note here that estimates of J\ and J 2 in 
terms of the downfolded ti values^ are within the same 
order of magnitude as those obtained here from the sus- 
ceptibility comparison. 
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FIG. 5: Temperature dependence of the susceptibility 
(emu/mol) at low temperatures. 

The partially frustrated model is invariant under trans- 
lations i — > i + 3 [i — 1, . . . , 9) and corresponds there- 
fore in the low-energy sector to a three-site ring with a 
four-fold degenerate ground state, i.e. it is a chirality- 
degenerate spin-doubled. The susceptibility of the par- 
tially frustrated model therefore diverges at low temper- 
atures as 1/T (illustrated in FIG. 0). Below ~ 10 K the 
experimental susceptibility increases slower than the cal- 
culated one suggesting the absence of a spin gap. This is 
the range of energy where the inter-ring couplings, which 
have been neglected so far in our discussion, will start to 
be important ( our ab initio t± — 0.03eV corresponds to 
a Jj_ ~ 10 K) 

We consider in the following the contribution of the 
inter-ring couplings. Our ab initio results show that - 
though small- there are various contributions, ti±, t 2 ±, 
t' 1± and t' 2± , to the inter-ring coupling (see TABLEQJand 
FIG. [21 which will be important at low temperatures and 
can compete against each other. These hoppings can be 



taken as a reference for the important exchange integrals 
to be considered, i.e. J^j_ and J' i± . Assuming antiferro- 
magnetic couplings, consideration of only n.n. Jx will 
generally induce a spin-gap^ which will scale with the 
inter-ring exchange coupling. On the other hand, the ex- 
istence of two competing exchange integrals Jj_ and J' ± 
which induce inter-ring frustration changes this scenario 
completely. As argued previously, the partially frustrated 
nine-site ring model corresponds in the low-energy sector 
to a three-site ring. We can therefore consider at low 
temperatures a model of three-site rings with frustrated 
inter-ring couplings which maps the main features ex- 
tracted from our ab initio calculations. The recent den- 
sity matrix renormalization group study by Liischer et 
al& on a tube of 3-site rings with frustrated inter-ring 
couplings showed that the spin gap diminishes if inter- 
ring frustration is present. For perfect frustration, i.e. 
J^=J'j_ the spin gap goes exactly to zero since then we 
are left with an effective 1-dimensional S=l/2 Heisen- 
berg model. Our ab initio hoppings t± and t' ± are of the 
same strength what strongly points to the condition of 
almost complete frustration of the exchange integrals and 
therefore to a spin-gapless behavior as has been observed 
exp er iment ally&. 

Finally, it is worth mentioning here that the partially 
frustrated nine-site ring model, is of considerable interest 
by itself. For a' = 1 and large number of sites it is 
geometrically equivalent to the spin-1/2 diamond chain, 
which corresponds in the low-energy sector to a mixed 
spin-1 spin-1/2 chain2£*2&. Our result a' = 0.7 indicates 
that Na2Va07 is close to the spin-diamond geometry. 

We acknowledge useful discussions with V.N. 
Muthukumar and J. Richter and the support of 
the Deutsche Forschungs Gemeinschaft. We thank the 
MPI-India partner group program for collaboration 
support. 



1 see for a review P. Lemmens, G. Giintherodt and C. Gros, 
Physics Reports 375 1 (2003). 

2 P. Millet, J. Y. Henry, F. Mila and J. Galy, J. Solid State 
Chem., 147 676 (1999). 

3 M.-H. Whangbo, H.-J. Koo, Solid State Comm. 115, 675 
(2000). 

4 V. Subrahmanyam, Phys. Rev B 50, 16109 (1994). 

5 H.J. Schulz in "Correlated Fermions and Transport of 
Mesoscopic System" Ed. T. Martin, G. Montambaux, J. 
Tran Than Van ( Frontiers, 1996). 

6 K. Kawano and M. Takahashi, J. Phys. Soc. Jap. 66, 4001 
(1997). 

7 H-T Wang, Phys. Rev B, 64, 174410 (2001). 

8 J. L. Gavilano et. al., Phys. Rev. Lett. 90, 167202 (2003). 

9 A. Liischer, R. M. Noack, G. Misguich, V. N. Kotov, and 
F. Mila, Phys. Rev. B 70, 060405 (2004). 

10 O.K. Andersen and T. Saha-Dasgupta, Phys. Rev. B62, 
R16219 (2000) and references there in. 

11 JP. Perdew, S. Burke, and M. Ernzerhof, Phys. Rev. Lett. 



77, 3865 (1996). 

O. K. Andersen, Phys. Rev. B, 12, 3060 (1975). 
P. Blaha, K. Schwarz, G. K. H. Madsen, D. Kvasnicka and 
J. Luitz, WIEN2k, 2001. ISBN 3-9501031-1-2 
We checked convergence of the band-structure results ob- 
tained with the TB-LMTO scheme -where the downfolding 
procedure is implemented- versus the results from the full 
potential LAPW scheme. 

We note that though the GGA calculations predict a 
metallic behavior for this system, since we have half-filled 
V—xy bands (corresponding to a V 4+ in 3d 1 configura- 
tion), it is to be expected that explicit inclusion of correla- 
tion effects through, e.g. an LDA+U calculation, will open 
a gap at the Fermi level and define the system as a Mott 
insulator. 

See references in O. K. Andersen, T. Saha-Dasgupta and 
S. Ezhov, Bull, of Mater. Sci. 26 19 (2003). 
The first neighbor shell of intra-ring V-V pairs are of dis- 
tances: 2.91 A, 2.98 A and 3.06 A. The next set of intra- 



■5 



ring interactions are of distances 3.71 and 3.55 A. 

Y. Ueda, Chem. Mater. 10, 2653 (1998). 

R. Valentf and T. Saha-Dasgupta, Phys. Rev B, 65, 

144445 (2002). 

R. Valentf, T. Saha-Dasgupta, J.V. Alvarez, K. Pozgajcic, 
C. Gros, Phys. Rev. Lett. 86, 5381, (2001). 
V. V. Mazurcnko et. al, Phys. Rev. B 66, 081104(R) (2002) 
and references therein. 

H. Rosner et. al, Phys. Rev. Lett. 88, 186405 (2002), and 
references therein. 

To be precise, one should add to J AF the J FM contribu- 
tion as well, which may be estimated by evaluating the 
Coulomb exchange integrals between two Wannier func- 
tions centered at two neighboring sites. However, since we 



intend to extract from the ab initio calculation only the in- 
formation of the important interaction pathes which should 
define the relevant exchange pathes, we have not included 
the FM contributions in our initial guess. 
For a U= 4 eV typical for vanadates and neglecting V 'inter, 
we obtain effective Ji= 285 K and J^= 211 K. The ferro- 
magnetic contribution -not included here- will renormalize 
this estimate. We note that the experimental Ji is about 
200 K. 

N. Fukushima, A. Honecker, S. Wessel and W. Brenig, 
Phys. Rev. B 69, 174430 (2004). 

S. K. Pati, S. Ramasesha and D. Sen, Phys. Rev. B 55, 
8894 (1997). 



